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Abstract 

We comment on the derivation of the main equation in the bounded 
confidence model of opinion dynamics. In the original work, the equation 
is derived using an ad-hoc counting method. We point that the original 
derivation does contain some small mistake. The mistake does not have 
a large qualitative impact, but it reveals the danger of the ad-hoc count- 
ing method. We show how a more systematic approach, which we call 
micro to macro, can avoid such mistakes, without adding any significant 
complexity. 
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1 Introduction 

One of the most famous mathematical models of opinion dynamics is presented 
in [5], in which the agents have continuous opinions and adjust them based on 
random encounters with other agents given that their opinions are below a fixed 
threshold d. Using a Taylor approximation and an ad-hoc counting method, the 
authors arc able to derive an equation for the evolution of the distribution of 
opinions for small values of d. 

We show that the original derivation does contain some small mistake. The 
mistake does not affect severely the model, as the qualitative behavior remains 
unchanged, but some values like the speed of convergence or the n-th order 
moments are impacted. This reveals the danger of the ad-hoc counting method. 

We then show how the equation can be derived in a safer way, using a more 
systematic approach, which we call "micro to macro" . It consists in explicitly 
modelling the microscopic system before going to the fiuid limit. The key tech- 
nical tool is the drift equation. We illustrate that this is not more complex than 
the ad-hoc approach, and can be more easily validated. 



*J. Gomez-Serrano's research was done while being an exchange student at EPFL. 
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2 The Model 



Wc consider a population of N peers, each of them having an opinion Xi G [0, 1] 
about some common subject. At every time step k € N, two peers are selected 
uniformly at random. If the distance between their opinions is sufficiently small, 
each of the opinions is modified, whereas if the opinions are far away they don't 
change. This is controlled with the uncertainty parameter d. Mathematically, 
we have that: 

• If \xi — Xj \ > d the opinions remain unchanged. 

• If \xi ~ Xj \ < d the opinions change the following way: 



Xi{k + 1) - Xi{k) = fi ■ {xj{k) - Xi{k)) 
Xj{k + 1) — Xj{k) = /i • {xi{k) — Xj{k)) 



where is a parameter (adaptation capacity) between and i. 

3 The Ad-Hoc Counting Method in P [T0| 

The final result is an equation for pit, x) defined as the PDF of opinions at time 
t: 

-^-h+h (1) 

where Ii and I2 are described next. In [S] only the final result is stated without 
proof; the proof is available in [TD] , but in French only. Therefore, in the rest 
of this section we give a literal translation into English of the relevant parts of 

m- 

A simple enough calculation shows that, at the limit of the small values of 
d, the system has the tendency to amplify the irregularities of the population 
density according to their opinions, p{x). During an elementary timestep cor- 
responding to one interaction, the variation of p{x) can he seen as the sum of 
two contributions: 

A negative contribution corresponding to the probability for an agent of opin- 
ion X to interact and modify his opinion: 

h = -p(x) I p{x)p{x + y)dy (2) 

J-d 

A positive contribution I2 corresponding to the probability that an agent of 
an initially different opinion has opinion x after the interaction. 

For this term, the obtaining of the normalizing constant is less trivial. The 
solution I have used is to consider a model in which the displacement of x -\- y 
while interacting with x -\- z is done around p{y — z) with a gaussian probabil- 
ity. Letting the length of the gaussian go to zero, one gets our initial model 
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conveniently normalized. The positive term is therefore calculated by the double 
integral: 



h =+ J dy ■ p{x + y) j dz ■ p{x + z)5{{x + y)+ p{{x + z) - {x + y)) = x) 

= + j dyp{x + y) J dz ■ p{x + z)S{y + p{z - y)) 

= + [ dy-p{x + y) [ dz ■ p{x + z)—S(z + -y) 

J J M 



- p{x + y) ■ p{x H y)dy 



(3) 



For small d, one can replace in the integrals the terms in p by their Taylor 
series around x (limiting to second order in y): 



Sp^ - J ^p- (^p + yp + y P"^ dy 



fj,d 



-I ip + yp' + \p' 



\p' + I ( ] P" 1 dy (4) 



P 



P 



where all the terms in p are evaluated at x and ' represents the partial derivation 
with respect to x. Continuing the calculation, one obtains: 



--{2pd)p^ + - 



P 



P 



5p=-2dp'-^p" 



2-^ (5) 



After simplification, and supposing that the characteristic interaction time 
is T, one finally finds at the first non-zero order: 



dp d^ d^P^) 
- = -.p.{p-l).-^ 

There are three errors that alter the final result published in [5]. 



(6) 



(i) Equation ^ contains an extra term p{x) and equation ([5]) lacks a factor 
p in the second term. 

(ii) In both Ii and I2 there is a multiplying factor 2 missing. Correcting both 
mistakes Ii should read: 



h = -2p{x) I p{x + y)dy 
J-d 



(7) 
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and l2'. 



h = - I p{x + y)p (x + ^ — ^-y] dy (8) 

The error here is a division by 3/2 instead of 2. Note that the error 
is propagated to [H] and, to the best of our knowledge, has not been 
corrected anywhere else. 

It is clear that the first correction is a typo, as it is corrected afterwards 
and the error is not propagated. Regarding correction (iii), it is a calculation 
mistake. Although the system's qualitative behavior isn't affected by this one, 
it might change the convergence speed of the system. Correction (ii) affects the 
time scale of the model by a factor of 20. This is a common difficulty in the 
ad- hoc counting approach, which can be circumvented if, as we discuss next, we 
specify the microscopic model in detail. 

The equation should read then, following the notation from [TU] : 

'^^ - ' p[ ^ ^,t]piy,t)dy^2p{x,t) p{y,t)dy. 



PJx — fid \ P / J x — d 

(10) 

or, in a more compact (symmetric) way, as in a Boltzmann-like equation 
ll],for/i^i: 



^ = ^1.(2.-.) \ 2,-1 'Va 2,-1 

^2p{x,t) / p{y,t)dy (11) 

J x — d 

Note that this mistake has also propagated to [T], in which it affects the 
calculation of the evolution of the second order moments, as it should read 
M2{t) = M2(0)e~*^°* (there is a factor 2 in the exponent missing), and from 
there to several papers such as [3] or [5]. However, the derivation is also cited 
in [9], where it is partly corrected. 

^The system is described by "A chaque pas de temps, deux agents sont tires au hasard 
dans I'ensemble de la population" (At every timestep, two agents are randomly chosen from 
the population), which differs from what is computed in the sequel: the variation of p{x) 
should be calculated as the positive and negative contribution (which are correctly written) 
for the case when the peer with opinion x is the first element of the randomly chosen pair, 
and again, the same contribution (because of the symmetry of the problem) for the case when 
he is the second element of the randomly chosen pair. 
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4 An Alternative Method Based on Micro-to- 



Macro 

In this section wc describe how the evohition equation can be obtained in a 
safer way, using a micro to macro approach. The idea [5] is to first write the 
"drift equation" for the microscopic system, i.e. for the discrete time system 
defined in Section [5] where the number of peers is finite. Then, in a second step, 
the evolution equation for densities as in Eq. (jlOp is obtained in the hmit of 
infinitely many peers. 

In the microscopic system, the density is not a proper one, but a mixture 
of Dirac masses. Instead of manipulating densities, it is more convenient (and 
faster) to manipulate the integral of arbitrary test functions against the density. 
Therefore, we proceed as follows. 



Step 1 (Micro) Let {t) be the occupancy measure of the microscopic sys- 
tem with N peers, defined as M^{t) — J2n=i ^x'^{t)- Let h be any bounded, 
measurable function defined on the space of one peer, i.e. on [0, 1]. The "drift" 
of the microscopic system is the operator defined by 

g^(/i)(i.) =E{h{M'' (t + 1)) -/i(A/^(t))| A/^(t) = 

where the starting condition must be of the form i/ — ■ With 

some elementary manipulations we obtain: 

N(N - 1) TV ^ iH^J'Xrn + (1 - fJ-)x„) + h{fiXn + (1 - tJ.)Xm) " HXm) ~ h{Xn))l{\xr,-x^\<d} 

^ ^ m<n 

= j^2(^ ^ I) X! ('^(/^^'" + (1 ^ f-^)^n) + HjlXn + (1 - ^J)x,n) - h{x„i) - h{Xn))l{\x^-x^\<d} 

+ ^2|^jV - n ^ {HP-Xm + (1 - lj)Xn) + HfiXn + (1 - m)-^™) - h{x„t) - h{Xn))l{\ 
^ ' n<m 

= j^.2(^ _ I) ^{HfJ'Xm + (1 - fJ-)Xn) + KjXXn + (1 " y)Xm) " Kx^) - h{Xn))l{\x„^x^\<d} 

{h{fiy + (1 - fi)x) + h{nx + (1 - Ai)y) - h{x) - h{y))l{\^_y\<j^~^dv{x)dv{y) 

[0.1]2 

{h{^ix + (1 - ^i)y) - h{x))l{\^_y\<a}dv{x)dv{y) (12) 



1 



TV- 1 
2 

TV- 1 



[0,1]2 



Step 2 (Micro to Macro) Assume that when the number of peers TV grows 
towards oo the occupancy measure has a limit, say vt- We must also re- 
scale time by a factor of N , as we see next, so that [k) w Vk/N- Whether 
this holds, and why, is an entirely different problem, which is not addressed in 
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[SI [TU] and which do not address either her^. Our purpose in this letter is solely 
to obtain a safe method for writing the equation. The method in e.g [3] says 
that, if we re-scale time by 1/iV, then the limit must satisfy 



d{h,vt) 
dt 



= lim Ng''{h){vt) 



if the limit exists. Intuitively, this is because in the limit we have {h,i>^^j_) 
{h,vt) ~ Q^{h){vt)- By Eq.dEl), we obtain: 

d{h,vt) _ 



dt 



[0,1] 



{h{^ix + (1 - ^i)y) - h{x))l{\.j,_y\<a^}dv{x)dv{y) (13) 



which is the main equation. To derive an equation for the density, we write 
vtidx) = p{t,x)dx, vt{dy) — p(t,y)dy and re-arrange Ea. ([T51) to obtain 



[04] 



hix)^dx 



[0,1] 



-pix,t) 



h{x) 



1 {x-{l-p)y 



,t)p{y,t)dy 



x+d 



Piy,t) dy 



x — d 



dx 



Eg. (fT2l) follows by identification of ^^fc^-* to the term between square brackets. 



5 Conclusion 

We pointed out that there were some mistakes in the derivation of Equation 
([6]) in [5]. The mistakes do not affect the qualitative behavior of the system 
and were partly corrected in later works, however, their existence suggests that 
the ad-hoc counting method for deriving such equations is error prone. We 
described how a more systematic method can be used to avoid such pitfalls. 
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